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ABSTRACT 


A survey of the literature on pressure, temperature and 
foreign gas inlet distortions is made. For interpretation of 
the influence of inlet distortion on engine stability, a break 
is made with distortion methods which use temperature or pres- 
sure distortion maps. The distortion in the form of tempera- 
ture contours is transformed to vorticity using appropriate 
equations derived for both incompressible and compressible 
flow. Crocco's theorem has been expanded in scope to include 
multispecie, nonuniform, perfect gas mixtures. The derivation 
of the equations and specialization to cases of interest in- 
volved a thorough order of magnitude analysis. A numerical 
solution of Poisson's equation as formulated by Squire and 
Winter to describe secondary flow arising from convected vor- 
ticity is presented. Deficient areas of knowledge relating to 
the general temperature/foreign gas problem are outlined. A 
conceptual model involving compressor stage transfer functions 
and feedback mechanisms is introduced and discussed at several 
points. Although not essential to understanding the vorticit} 


stall logic chain, this model is a useful conceptual tool and 


may be worthy of additional study. 
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SYMBOLS 
Cy specific heat for a mixture of perfect gases 
Cy Specific heat for a single perfect gas 
D matrix of finite difference coefficients 


FG foreign gas - variations of Y and R 
G forward transfer function 

H fecdback transfer Tunction 

h specific enthalpy 
M 


Mach number 


m mass/unit volume 
ig) normal direction in natural coordinates 
; P pressure 
7 q scalar velocity 
r radial displacement in cylindrical coordinates 
Ry universal gas constant 
R R,/m . 
s specific entropy 
T temperature 


TFG temperature and foreign gas - variations in T, yY and R 


U velocity 


Vv velocity perturbation 
° W velocity perturbation 
; y coordinate direction perpendicular to cascade span 
1 


There is a variation of this notation in Appendices” 
A,B,E, and G. Where such occurs, it is clearly indicated. 
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vA axial displacement in cylindrical coordinates 

r circulation 

Y ratio of specific heats 

A finite difference 

€ flow turning angle Eq. (19), order of magnitude 

8 circumferential displacement in cylindrical coordinates 

E streamwise vorticity 

p density 

co) angle between principal normal and the normal to the 
Bernoulli surface 

yp stream function 
total flow vorticity 

W component of flow vorticity 

Subscript 

r radial component 

t stagnation quantity 

Zz axial (streamwise) component 

8 circumferential component 

1 inlet 

2 outlet 


Superscript 


t 


mean value 


dimensionless quantity 
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I. INTRODUCTION 


Turbojet stability is often defined in an inverse sense - 
in terms of Reo tbal Ly. More specifically the turbojet 
stability is defined by the components most susceptible to 
instabilities - fans and compressors [Ref.1]. For a compres-~ 
sor characteristic (pressure-ratio versus corrected flow at 
constant speed) the border between stable and unstable opera- 
tion is the point of maximum pressure ratio. At this point 
the blade loading is such that the blades are aerodynamical- 
ly stalled. The resulting unsteady phenomenon is called 
surge. The difference between the pressure-ratio at the 
point of surge and the pressure-ratio at which the compres- 
sor is operating is defined as stability margin. 

A common factor causing compressor surge is non-uniform 
inlet flow. The non-uniformities at the inlet are popular- 
ly termed inlet distortions. 

The stability margin of current high Seeetfic energy 
gas turbine engines is at what might be termed a practical, 
if not critical, minimum. A result is that engine inlet 
distortions can readily produce instability in the fan 
and/or compressor. For certification purposes, the instabil- 
ity problem can be avoided by increasing stall margin 
which generally is at the expense of ‘engine performance. 

The most prevalent inlet distortions causing instability 


are total pressure and total temperature. These can be 


lal 





further categorized [Ref. 1] as: 


Inlet Pressure Spatial Distortion 
Inlet Pressure Fluctuations 

Inlet Pressure Ramps 

Inlet Temperature Spatial Distortion 


Inlet Temperature Ramps 


Pressure distortions are seemingly the most common in 
occurrence and thus have been the object of widest study. 
References 2-14 provide insight and additional literature 
sources for further review of this problem. 

Temperature distortions in turbojets received some early 
investigation [Refs. 15, 16, and 17] which centered on the 
transient category. More recently Rudey and Antl [Ref. 18] 
investigated both spatial and transient temperature effects 
on a turbofan engine compressor system. 

Other recent investigations have concentrated on the 
inherently related areas of steam and armament exhaust in- 
gestion. Tomassetti [Ref. 19] performed extensive work 
on the TF30 turbofan engine. He concluded that oie primary 


factors causing instability were: 


Inlet Temperature Distortion 

Non-homogeneity of Inlet Gas Properties 

Ingestion of Liquid Condensate 
Tomassetti developed a steam Envirommental Factor as means 
of comparing results. Reference 20 modified the Environ- 
mental Factor to accommodate additional variables. Also 


investigating the TF30 steam problem were Mallett and 
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Parcells [Ref. 21] and Worobei [Ref. 22]. They reached 
conclusions similar to Rudey and Bab Ls Additionally, 
Mallett and Parcells noted a dependency of stability margin 
on the individual compressor design. 

Armament gas ingestion has been the subject of several 
studies at the Naval Air fopainien Test Center (NAPTC). 
Rich and Reale [Ref. 23] investigated the effects of sea 
level missile exhaust gas ingestion on the J52-P-6A engine. 
The missile exhaust gas simulator used in NAPTC testing is 
described by Rich (Ref. 24]. Reference 25 publishes the 
results of rocket gas ingestion by the TF34 engine, which is 
used on the S-3 Viking and A-10A attack aircraft. af 

To gain further insight into the temperature problem, 
this paper will develop a mathematical model to describe 
the inlet distortion in terms of vorticity. This is a new 
approach in contrast to the classical approach of a total 
temperature map. It is an aid to a better examination of 
the flow mechanisms involved with inlet temperature 


distortions. 
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II. TEMPERATURE DISTORTION 


A. TEMPERATURE SOURCES 

Engine instability due to temperature or foreign gas 
distortion has not been widely studied relative to distor- 
tion from variation in stagnation pressure. Perhaps the 
reason is the complexity of the investigation due to the 
type of temperature distortion sources: i.e. catapult 
steam, Fig. 1; rocket or missile armament exhaust gas, 
Fig. 2; spent gun propellants; exhaust of other aircraft; 
and uneven temperature due to heating of air above runways 
and taxiways. The steam ingestion problem is of importance 
to carrier launched aircraft, especially those with low- 
mounted inlets on the aircraft centerline, 1.e. A-7. The 
ingestion of armament gas is mainly encountered by tactical 
aircraft. In addition to the detrimental effects of engine 
stall caused by the swallowing of rocket plumes, explosions 
may occur in the engine (the exhausts of rocket plumes are 
usually fuel rich). In formation flying or during aerial 
refueling an aircraft may be subjected to the exhaust gases 
of the leading aircraft. Uneven heating of air can occur 
at airports; air above runways, particularly blacktop run-~ 
ways, can be considerably hotter than the ambient air. 
During taxi or takeoff the aircraft engine can swallow 
streamtubes of air with this uneven ‘temperature distribution. 


It has been found that some aircraft with ejector nozzles 
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Figure 1. Catapult Steam Ingestion 


(adapted from Ref. 21) 
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Figure 2. Missile Exhaust Gas Ingestion 


(adapted from Ref. 24) 
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may experience temperature-induced stall while taxiing down- 
wind. Wind blowing up the tailpipe carries air along the 
secondary air passages where it is heated by hot engine 
components and thence into the compressor. 

Granted, sources of temperature distortion appear in 
limited situations. This does not mean the problem is in- 
Significant. Certainly, an engine stall during an aircraft 
carrier catapult launch is critical. Likewise engine in- 
stability during air-to-air or air-to-ground combat maneu- 


vers for aircraft is serious. 


B. CHARACTERISTICS OF TEMPERATURE SOURCES 

From the discussion of possible temperature distortion 
sources in the previous section, there are three character- 
istics that should be noted: 1. The distortion source is 
independent of the affected inlet-engine system. 2. The 
distortions may be both spatial and/or transient. and 
3. Foreign gases may also be a distortion source. A brief 
discussion of each of these characteristics is warranted. 

Necessarily, some definitions are required prior to the 
discussion. When the word "temperature" is used, it refers 
to total temperature. If static temperature is being dis- 
cussed, it will be clearly indicated as such. The phrases 
"temperature and foreign gas" and "foreign gas" will be 
used with sufficient frequency to warrant the abbreviation 


TFG and FG respectively. A symbols sheet on page 8 provides 


a complete listing of definitions. 





The significance of the independence of the source lies 
in the possibility of eliminating the distortion instabili- 
tiles by removal of the source. A major part of the effort 
in solving the A-7 steam ingestion problem involved minimiz- 
ing the source of steam [Refs. 19 and 22]. The aircraft 
exhaust problem can be minimized by judicious pilot tech- 
niques. Conversely, the removal of the source of armament 
gas is unlikely to provide a satisfactory solution. In 
examining the TFG problem, independence of source and engine 
Should be kept in mind as a potentially controllable variable. 

When examined from the P. or Ty map viewpoint, spatial 
distortions can be divided into circumferential and radial 
distortions, as illustrated in Figs. 3a and 3b. In the fan 
engine, the circumferential distortion causes a shift in 
the surge line of the fan and compressor [Ref. 1]. Figure 4 
is representative of this shift. Radial distortion can 
have similar effects except that the fan can cause attenua- 


tion of the distortion leaving the fan resulting in a lesser 





Circumferential Radial 


Distortion Distortion 
(a) (b) 


Figure 3. Spatial Distortions 
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Figure 4. Effect of Spatial Distortion 


(adapted from Ref. 1) 
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distortion level [Ref. 1]. In practice the distortions 

are composed of elements of both [Ref. dal ce 
x Transient or ramp temperature distortions may result 

in compressor response as shown in Fig. 5. Analysis of this 

problem involves unsteady aerodynamics and is discussed in 

a later section of this thesis. 

The effect of the presence of FG is mainly in changing 

_— thermal and caloric properties of the mixture .> These 
e gas properties, ratio of specific heats (y) and gas constant 
(R,/m) then become variables of interest. Any model of TFG 
distortion (e.g. due to steam or rocket exhaust ) must account 
for changes of these properties. The instance of no FG 


(e.g. hot runway) then becomes a special case. 


LOCUS OF INSTANTANOUS 
OPERATING POINTS AFTER 
TEMPERATURE CHANGE 







«OPERATING LINE 








PRESSURE | CONSTANT SPEED——- 


LINE LOSS IN SURGE MARGIN 


CORRECTED AIRFLOW 


Figure 5. Effect of Temperature Ramp 


(adapted from Ref. 1) 


1additional effects in the form of diffusion, condensa- 


tion and chemical reactions can be present. These phenomena 
will not be considered by this thesis. 
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Lil. JNU DISTORTION AND VORTICITY - 


The efforts to solve the distortion-induced insta- 
bility problem have in the past been studied independ- 
ently by one of two design groups - inlev or 
blade-cascade. The inlet group concentrates on analy- 
sis of instabilities utilizing events occurring at the 
compressor face i.e. inlet distortion. Their interpre= 
tation must be based on the distortion's influence on 
compressor internal aerodynamics. The blade-cascade 
group analyzes blade and cascade response to various 
flow situations. In terms of solutions to current 
distortion-instability problems, the inlet group 
efforts are vital. In achieving the ultimate goal of 
accurate prediction of flow of an arbitrary compressor 
with arbitrary inlet flow, it seems the blade-cascade 
group efforts offer most promise. A germane quegt ton 
is "Can the efforts of both groups be interfaced?" 

The approach of the inlet group to analyzing com- 
pressor inlet distortions involves the-use of a com-— 
pressor map (a series of maps for time-dependent 
situations) of inlet isotherms or isobars. An inlet 
map of the second stage would show the inlet distor- 
tion as modified by the aerodynamic influence of 


the first stage. Mokelke [Ref. 26] found through 
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experiment that the inlet guide vanes had very little... 

effect on distortion whereas noticeable alteration was pres- 
ent after the first stage. The term aerodynamic influence 
encompasses the phenomena of shed vortices due to non-uniform 
and/or unsteady lift and cascade secondary flows. This 
aerodynamic influence is sveeieeis the object of study of 

the blade-cascade group. It actually can be viewed as a 
transfer function. If one knew or could predict the trans- 
fer function of each stage, a large step would have been 
taken toward the previously mentioned goal. To begin the 
a function process, inlet distortion maps, temperature 
and/or pressure, must be converted into the conceptually 
compatible aerodynamic influence quantities. 

Imagine a compressor stage preceding the first stage of 
an actual compressor. For identification call it the zeroth 
stage. The function of this imaginary stage is to transform 
an initially uniform flow into a flow with aerodynamic dis- 
tortion patterns identical to the actual TFG (or pressure) 
distortion pattern. In essence the first stage of the com- 
pressor is affected by the aerodynamic influence of the 
zeroth stage. Thus the zeroth stage offers a way to view 
inlet distortion in the same’ physical sense it affects the 
first stage. It is also grounds for considering an interface 
between the inlet and blade-cascade group. The fluid mecha- 
nism involved is vorticity. 

The logic chain into which vorticity fits involves: 


1. conversion of distortion maps into vorticity maps; 
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2. evaluation of vorticity maps to determine changes in 
compressor blade loading; 3. utilization of blade loading 

to evaluate stall criteria; and 4. decision as to whether 
stall will or will not occur. The notions of aerodynamic 
transfer function and zeroth stage are not essential in under- 
standing the logic chain. They are introduced as convenient 
‘ways of visualizing distortions in terms of vorticity. In 
fact, the zeroth stage is analogous to Marble's actuator 

disk model [Refs. 27 and 28] which describes the flow through 
a single blade row. 

Before developing the partiality relationships, a discus- 
sion of the elements of aerodynamic influence is considered 
necessary. The concept of circulation about a wing is well- 
known. If the lift changes, the circulation changes, result- 
ing in shed vortices. The unsteady and cascade secondary 


flow effects are considerably more complex. 


A. BLADE AND CASCADE UNSTEADY FLOW 

Even if the flow entering the first stage is steady, it is 
not difficult to imagine a circumferential distortion causing 
unsteady flow over a rotor blade. In fact, Carta [Ref. 29] 
found that if the blade-passing reduced frequency was greater 
than 0.1, the flow is considered to be unsteady. During steady 
operation of a compressor the interference effects of the 
cascades are unsteady [Ref. 30]. 

Sernamdetlick [Ref. 31] has divided the unsteady phenomena 
into two groups: 1. purely aerodynamic unsteady phenomena; 


and 2. phenomena requiring coupling between the aerodynamic 
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and mechanical characteristics (aeroelastic effects). To 
limit the scope of this discussion, the assumption of infi- 
nitely rigid blades is nade, 

It has been shown by Carta [Ref. 29] that unsteady flow 
over an airfoil delays stall which produces unsteady Jeet 
coefficients above those of steady flow (Figs. 6a and 6b). 
He also found that the rate of change of incidence angle 
had a major effect on the unsteady response (i.e. note the 
severity of dC / do in Fig. 6a). The overshoot of Cy is due 
to the large reduction of adverse pressure gradient [Ref.7]. 
Ericsson and Redding [Ref. 32] explained the overshoot and 
undershoot as being caused by: 1. the erréect of induced 
flow acceleration on the adverse pressure gradient; and 
2. the effect of the induced change of effective Reynolds 
number or turbulence level on the boundary layer. 

From basic airfoil theory, changes in lift require chan- 
ges in wing circulation. Therefore unsteady lift produces 
an unsteady vortex wake. These wortices influence the flow 
about the originating airfoil. Additionally, they are swept 
downstream where they contribute to the interference effects. 
The interference effects are not necessarily detrimental. 
Shorr and Reddy [Ref. 33] found that in comparison with a 
single airfoil the amplitude of the unsteady lift coefficient 


for cascades reduces more slowly with increasing frequency. 





leor a comprehensive treatment of aeroelastic vibration 
in turbomachine cascades, see Samoylovich [Ref. 31]. 
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Angle of Attack 
(b) 


Figure 6. Unsteady Normal Force Loops 


(adapted from Ref. 29) 
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They concluded that the trailing vortices of a cascade have 
a mutual cancellation effect. 

It should be pointed out that this discussion was based 
on unsteady effects due to pressure distortions. To readily 
generalize these effects to the TFG case may result in inac- 
curate conclusions. This is an area where experimental TFG 
data would provide an increased understanding of the pheno- 


menon of unsteady, TFG flow. 


B. CASCADE SECONDARY FLOW 

A definition of secondary flow in turbomachines is that 
it is the difference between the actual flow and the ideal- 
ized axisymmetric flow [Ref. 34].-. It has been estimated 
that for high specific energy compressors up to 80% of the 
loss in total pressure is attributable to secondary flow 
effects [Ref. 35]. Figure 7 is a qualitative representation 
of most of the components of secondary flow and the resulting 
vorticity. It is easy to see why accurate prediction of the 
aerodynamic influence of a cascade alone has not been 
achieved [Ref. 30]. In considering the effects of inlet TFG 
distortion on secondary flow, attention will be focused on 
those components of vorticity modified by turning the flow.+ 

Squire and Winter [Ref. 37] developed a model for steady, 
incompressible and inviscid flow which treated the secondary 


flow as a linear perturbation superposed on the main flow. 


leor a state of the art report on the analytical treat- 


ment of secondary flows, see Woods [Ref. 34]. 
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Figure 7. Secondary Flow and Vortices in an Axial Flow Compressor Rotor 


(adapted from Ref. 36) 








Applying linear theory, they were able to show the genera- 
tion of a streamwise vorticity component when inlet vortici- 
ty normal to the streamtube is turned by a passage. Haw- 
thorne [Ref. 38] extended the turning effects to inelude 
trailing filament vorticity and trailing shed vorticity. 
The former is due to stretching of the vortex filaments be- 
tween the upper and lower stagnation streamline (Fig. pe 
The latter is due to the change in circulation along the 
blade. The net value of secondary vorticity is the sum of 
these effects. Section VI will present an analysis of the 
Squire and Winter effect. 

Again, it is worthy of note that these secondary flow 
effects are due to vorticity. A search of the available 
literature indicates a need for study of secondary flows 


resulting from TFG distortions. 
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Figure 8. Vortex Filament Passage through a Cascade 


(adapted from Ref. 13) 
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- IV. VORTICITY FROM TFG DISTORTION 


A. DERIVATION OF MODEL 

By proper choice of coordinate system and some simplify- 
ing assumptions,it is a straightforward matter to develop a 
relationship among vorticity, temperature, and pressure 
which is valid for a mixture of perfect gases. Such a devel- 
opment is presented in Appendix A. For a dual perfect gas 


system the result in natural coordinates is 


P 
Uw = a “ : a - $e ‘ (1) 

- Since this equation does apply to a mixture, enthalpy is 
defined as h = oe where e is the specific heat of the mix- 
ture. In the majority of cases involving TFG distortion, 
the inlet is exposed to a non-homogeneous mixture (i.e. 
there are spatial gradients in bie components of the mixture). 
Thus there will be gradients in Coe This clearly differen- 
tiates the case of one component (i.e. air) and two compo- 
nents (i.e. air and steam). 

Equation (1) can be derived in a much more general for- 


mat. In vector notation the continuity equation is 
+ ¥+(pU) = 0. 


For inviscid flow the equation of motion is 


pu. a0 ~ 2 Pee 
pee et. WH) T= - FUP. (2a) 


ct 
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A vector identity is 


GU) = (UV) + Ux (Wx) 


where ¥xG « ¥ is the vorticity vector. Combining Eqs. 


and (2b) yields 


C+ 


3 


a 


~ §xd + 5-H = - > VP. 
From the definition of total enthalpy 

Vn, = ¥(20-0) + Vh. 
Introducing Eq. (2d) into (2c) results in 


Ux = Vh, $ : a. 2 tee ee, 


For steady flow this reduces to Eq. (1). Equation (3) 


the model to be used to determine flow vorticity. 


tions involved in deriving this model are: 


Inviseid Flow 
Perfect Gases 


Entropy Constant Along Streamlines. 


B. VORTICITY IN INCOMPRESSIBLE FLOW 


The isentropic relationship 


P 
t ss y-1 
pet yel M?) 


can be expanded to 


me 2 mM? Me 
H, SF = 5 p Ukl tp tay) oy 


which for small local Mach numbers reduces to 
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(2b) 


(2a) 


(2¢) 


(2a) 


is 


assump- 








“ 2 
ier. = %5pU~ . (4) 


Introducing Eq. (4) into Eq. (3) results in 


=} 
a aw te & _ 4-2 aU 
ui x = - VP, — Ww Vine + = (5) 


Note that the presence of FG leads to the gradient in densi- 
ty.+ For a single perfect gas in incompressible flow, the 
gradient in density disappears leaving vorticity as a func- 
tion of velocity, total pressure, and the unsteady effects. 
In terms of experimental data gathering, these factors seem 
well within the scope of present techniques. However, the 


addition of FG to the flow creates a complex relationship. 


C. VORTICITY IN COMPRESSIBLE FLOW 


Expanding Eq. (3) 


> 
> A. dU 
Vee on, Vo. + 6.9. + = WP = = ey Se 
G70. C¥ly * 5 P TVC,, a st 

This reduces to 
Fadl 3 @o - Gr) +2 au 
Fx 8= (n, - TV, + O(N - a +5 Ye+5- (6) 


For the special case of one gas, the first term on the 
right side would vanish. But for the general case of both 
temperature and foreign gas,Eq. (6) applies. Appendix C 
contains a simple application to provide an example as well 


as an order of magnitude and sign verification. 





1 appendix B contains an example of how this result might 
be applied. 


spl. 





V. VORTICITY MAPS 


A. ORDER OF MAGNITUDE ANALYSIS 
The vorticity equation for compressible flow was found 


to be 


Gx 8 = (r, - TV, + clr, - Fr) +i Wee. (7) 


Realizing that a search for compressor inlet mass distri- 
bution data would probably be futile, the special case of 
temperature distortion was considered. The temperature dis- 
tortion experiment conducted by Rudey and Antl [Ref. 18] 
was found to closely approximate inlet temperature-only dis- 
tortion.+* Thus the flow can be considered to behave as a 


single gas which reduces Eq. (7) to 


oly 


a x = a a st 
Ux 8 esCval, vr) + ae: + (8) 


The number of variables can be reduced by eliminating Vr. 
The procedure involves introducing the combined first and 


second laws of thermodynamics 


Tds = dh - = aP 


into Eq. (8) resulting in 


alee 


: Ux = oT, - TVs + (9) 


1 appendix D contains an analysis to this effect. 








This is Crocco's Theorem. From thermodynamics 


> ¢C > > 
= —2 bey ee 
Vs a vt. p VP. ; (10) 


Combining Eqs. (9) and (10) leads to 
> oar RT a 
jx Gece tr ~- Pte + EP Wlte. 
p ot T t 
t t 
Using the definition of Cy and the isentropic relationship 


between T, and T, Eq. (10) can be rearranged to 


2 aa 
U x 4) = ae —_+,— VO, + RT ve + a“ (il) 
al, “4 — M? t 


For this analysis, an order of magnitude (ce) is defined as 


From Fig. 9 the maximum value for Ty, is 590°R. The mean 
temperature is approximately 520°R. Thus O(e) = 0.13. For 
the complete map, -0.13 £ € & 40.13. 

Returning now to Eq. (11) and examining the coefficient 


of vr, 


2 1 + a M2 2 2 


2 2 2 
YRM | 1 | ~ YRM F _ Yel ye ¢ a mn] +} 
For Mach numbers less than about 0.5,it is apparent that all 
terms beyond the first are O(e*) or less and are therefore 


neglected. This leaves 


_ YRM2 2, RT 30 
Ux $= = v1, + P VP, oS. (12) 
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Figure 9. Inlet Temperature Map 


(adapted from Ref. 18) 
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Equation (12) can be written in terms of eylindrical coordi- 


nate system components 





aT aP 3U 
HEN gp RE ae 
U4, ~ UZ % 2 or - rE. or + 3¢ » (12a) 
oT 9P aU 
_ YrmM? +t RT ‘ 8 
Ui, ~ Und, = Sp Be + Pr Be" Ut i120) 
ae aT P U 
9 3 
_ YRM? t RT °'t Zz 
Ub, > Ug%, Pe ae P, oz 7 ae (12c) 


To non-dimensionalize the vorticity equations, reference 
quantities were selected at the engine inlet. Based on data 
from Rudey and Antl [Refs. 18, 39 and 40], the quantities 
(denoted by a superscript bar) used in this analysis are: 
axe se 520°R; 2. T, 510°R; 3. Pr, 1.5 ft; and 4. M, 0.4. 


These reference quantities are used in the following defini- 


tions 
= ta =- PDP Pp 
U Ula P. PEPE 
wy = wta/P t = t'r/a 
= rr = i 
T, Tpit. Zz, Z'r 
ps Pip p Ss pir 


The superscript prime signifies non-dimensional quantities. 


Rewriting Eq. (lea) in non-dimensional form 


Since 


M2 = 1+ O(e?) = 1 
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the non-dimensional form becomes 


(3 A te eke ee, (13a) 


~- ! — 
Uae, UW 2 Ee Y 


Equations (12b) and (12c) are similarly transformed into 





oT! oP! ou! 
ae | .! : 6 
t te ces t , ee en a aie —asce 
Pte U4, ry 30 : yPir! re) "i at” (13b) 
and 
» oT! , oP! au! 
Uses - Uj = 5 agr + YET at sr (13¢) 


An order of magnitude for eT i /or' can be established us- 
ing Fig. 9. At the 8 o'clock position, the value for AT{ /Ar' 
is 2.5 or O(e°). Likewise, oT /(r'98) can be evaluated at 
the 2 Seioek position. There, ATI /(r'A8) was calculated to 
be 1.7 or O(e°). An approximation of oti /oz' can be ob- 
tained using data from Ref. 18. It was found that the engine 
inlet experienced temperature changes of up to 8500°R/s. 

By non-dimensionalizing this rate and dividing by UZ, the 
order of ott /oz" was found to be e. This involves the as- 
sumption that DT‘ /Dt! = 0. 

The components of vorticity are the unknowns to be 
determined. They are assigned O(e°) subject to later con- 
firmation. Order of magnitude estimates for the remaining 
terms are generally based on physical considerations, which 
are now discussed. 

The main flow can be considered to be almost totally 
axial. Thus Ups is essentially M which is O(c). The presence 


of my is caused by streamtube divergence from the heating 
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paseesas— If the engine inlet is sufficiently far enough 
downstream that the streamtube divergence approaches equilib- 
rium (i.e. static pressure across a plane normal to the 
flow is approximately constant), then U), is O(e*}. Whie 1s 
mathematically shown in Appendix E. If it is assumed that 
there is small duct swirl, then by applying an argument simi- 
lar to that of Uj, Ug is G(e*), 

Evaluation of aU, /ot! term was possible through use of 
Fig. D-3s Using finite differences between lines of con- 
stant velocity at constant z,a value for AUS /At' was computed 
as 0.006. Since the region of heat addition should also be 
region of maximum AU, /At' , assignment of 0(e*) seems Season 
able. An evaluation of aut / ot" and 9U'/3t' was conducted in 
Appendix E. As a result, they are assigned OLe™)z 

The term oPi/ oz! was evaluated in Appendix F and found 
to be O(e€7). The terms api /(r' 96) and oP, /or! are assigned 
O0(e2) based on information from Ref. KO. 

A summary of the preceding discussion is presented in 
Table I. It is reiterated that these results are estimates. 
Additionally they are applicable only to the flow system 
under analysis. 

Equations (13a, b and ec) are rewritten with appropriate 


orders of magnitude above each term 


lingine inlet hub effects on US are neglected, since the 
map of temperature is at a station ahead of the hub. 
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Table I. Relative Orders of Magnitude 


Variable Comments 


w! Assumed 
r 
Assumed 
Assumed 


Streamtube divergence due to heating 


Confined by duct 


Streamtube divergence due to heating 
No swirl added in bellmouth 
Swirl size of € would cause engine 


to stall 


Main flow 
Fig. 9 
Fig. 9 
Fig. F-3 
Fig. F=-3 
Appendix E 
Appendix E 
Ref. 39 


Ref. 39 


Appendix F 














a : c= ge «4 e< 
e*e E € y Ole , Pi = au, 
Tot o U) Ws = oo = + TE st + str ? (13a) 
0 2x10 - e" = e a 
€ ere 2 OT! . 8Pf a2 
6 
Uso, - Ula = dor ae + YET" sg * oe? CES) 
and : 
2.0 2 J e° - ad 
ete® ere 2 oT! eB) 80 
Uptige= Uge, =o gar ser Ter * HET oe 


Dropping terms and products of terms of order e* in Bae. 


(13a) and (13b) results in 


1 
Sites = BO al : 
So an an To (14a) 
and 
rt 
M2 Ott 
1 Rothe. ee 
Use), * set 95 . (14b) 
Since U! = M,the final result becomes 
oT! 
M t 
’ = -_ — 
"9 2 or! (15) 
and 
oT! 
M t 
ee 
Yr er' 30 (16) 
Note that for M = 0.4 and eT /or' = -2.5, the corresponding 


value for W6 is 0.5 which is not O(e°) as previously assumed. 
What is important is the product of Ur, is the same order 


as the product of (Mm?/2)(aTt/ar') as seen by 


(470.5) = (.08)02.5) 


34 2.45 = 1,22 2-85 


ee © € e 
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Thus the estimates of the individual terms may vary slightly 
from the actual values. But the hope is that the variance 
is small. Additional smoothing can be expected from taking 
products as was seen in this case. 

Equations (15) and (16) were employed to convert the 
temperature map (Fig. 9) into circumferential and radial 
vorticity maps (Figs. 10 and 11). The streamwise vorticity 
Cw) ) map could also be developed thru application of the 
relationship div(®) = 0. This enables one to map all three 
components of vorticity from a two dimensional map; there 


is, however, an ambiguity in the algebraic sign of Woe 


B. DISCUSSION OF RESULTS 

Before proceeding, further background concerning this 
particular temperature distortion induced stall is necessary. 
Rudey and Antl [Ref. 18] found that the engine stall occurred 
in the compressor and propagated forward. They attributed 
this to the fact that the hottest and most rapidly changing 
temperatures were concentrated in the hub region of the fan. 
This is the portion of the flow which entered the compressor; 
whereas, the cooler portions:-of the flow went through the fan. 

The vorticity maps are valid for only one point in time. 
Nonetheless, they do indicate that significant aerodynamic 
influence can be expected from the temperature distortion. 
Negative wy, causes the blade circulation to be increased thus 
increasing blade Lowes Positive we causes the inner por- 


tion of the blade span to experience increased flow 


4O 








Ps Figure 10. Circumferential Vorticity 
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Figure 11. Radial Vorticity 


He 








velocities and visa versa. In both cases Wn and W6 cause 
vortex shedding which is probably unsteady. Additionally, 
associated with W5 is a non-uniform velocity distribution 
-which along with We causes cascade secondary flow effects. 
Hence the vorticity approach converts a distortion into aero- 
Gynamic terms. 

An interesting feature of Fig. 10 is the large ring of 
WE at midradii extending from approximately 1 o'clock clock- 
wise around to 12 o'clock. In his study of vorticity due 
to total pressure distortions, Farmer [Ref. 13] found rings 
of We of approximately the same extent. As will be shown in 
section VI, WE has very definite effects on a cascade. It 
can therefore be stated that a ring of We» Such as in Fig. 10, 
is a good indication of instability. 

The lack of time history and engine performance data 
makes further analysis of these maps relative to their role 
in the engine stall difficult. In addition to the previous 
time aspect, time history is of importance hen. Gompeeine 
the transport time of the distortion between the point of 
measurement and the point where it induces instability. For 
example, if the instability occurred in the tenth stage and 
the data gathering took place at the engine inlet, the dis- 
tortion map of interest probably occurred during the previous 


ten to fifteen milliseconds. Rudey and Antl [Ref. 18] found 


the flow breakdown occurred approximately thirty milliseconds 


lRor a detailed analysis on inlet vorticity effects, see 
Farmer [Ref. 13]. 
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after the inlet temperature rise was detected. Figure 9 

was recorded about five milliseconds before stall. Since 
their inlet recording station was about a fan diameter 

ahead of the compressor, it would probably take fifteen to 
twenty milliseconds for an inlet distortion to reach the 

aft stages of the compressor. Based on this estimate, the 
inlet distortion causing instability occurred ten to fifteen 
milliseconds before that in Fig. 9. Therefore it is not 
known whether Fig. 9 is totally representative of the distor- 
tion map which was the advent of stall. According to the 
analysis in Appendix F one would anticipate little change 

in this time frame. 

Data poneeedine engine performance is necessary to deter- 
mine such factors as blade loading and estimates of flow 
turning angles. The former is required to predict effects 
on changes in blade circulation. It will be shown in the 
next section that the flow turning angle can be used at 


least qualitatively to estimate secondary flow losses and, 


more importantly, onset of stall. 





VI. SECONDARY FLOWS DUE TO INLET VORTICITY 


A. SECONDARY FLOW MODELS 
‘ The analysis of secondary flows has resulted in four 
flow models - inviscid, viscous, intermediate, and tip 
clearance [Ref. 34). The inviscid model is a considerable 
simplification. Even further sinod Ure wean are necessary 
if Hawthorne's [Ref. 41] result 


E E) 5 ; a ; de 
fe) _ (2) «= 2 [tre 21 ome an 


is to be solved. Hawthorne and Novak [Ref. 42] discuss | 
three classes of approximations required to produce solutions 
of Eq. (17). They suggested using the assumption of a pri- 
mary, irrotational flow which convects the vortex filaments. 
Thus, secondary flow is induced by the éonvected vorticity. 


The following discussion utilizes the suggested assumption. 


B. POISSON'S EQUATION 
Consider the system represented by Fig.12. From the 


definition of vorticity, the streamwise component is written 


_ ow ov 
eas Se (18) 
» For the assumption of steady vorticity superposed on the 
mainstream, the induced velocities are defined as v = a 
, € ow ov OW _ ‘: ; 
. and w= - dy such that continuity & + ay ois satisfied. 
y is the stream function. Combining these definitions re- 
. sults in Poisson's equation V*y = - & . The vorticity 
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Pigure 12. Model of Channel for Cascade Secondary Flow 
. (adapted from Refs 37) 
(e, ) is Soesdned by solving Eq. (17). The solution of 
Poisson's equation allows determination of the change in 


blade circulation ar. by [Ref. 42] 
= V5 a 
At. = , vdy . 


The change in lift can then be computed by 


AL = AT pu 
sm 


where Ue is the vector mean of the velocities of entry and 
exit, U1 and U2 [Ref. 43]. Through application cf momentum 
concepts, the change can be expressed in terms of kinetic 
energy. However, for the inviscid model the values calcu- 
lated would be qualitative in nature fRef. 44]. Thus a solu- 
tion of Poisson's equation only provides insight into 


secondary flow losses. 
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C. APPLICATION OF POISSON'S EQUATION 

In punee. to illustrate the effects of secondary flow, 
a solution to Poisson's equation was obtained. In so doing, 
it was necessary to determine values for Bo One approxi- 
mation often used in the literature is the Squire and Winter 


[Ref. 37] expression 


_ au | 
fy = — 2E ae ‘: (19) 


In Eq. (19), ¢ is the turning angle and U is the inlet 


velocity. The assumptions underlying this expression are: 


Steady, incompressible, inviscid flow 


Rectangular channel with 2 ee 


ax x 


Inlet velocity variation (vorticity) a function of z only 
Radius of curvature of the channel is large 


Induced velocities are perturbations on the main stream 


Despite these assumptions the results obtained by Squire 
and Winter were very satisfactory. In fact Eq. (19) can 
be applied to impulse turbines with reasonable accuracy 
[Ref. 42]. 

Due to the uncertainty in the results of Section V 
another flow situation was chosen. Farmer [Ref. 13] obtained 
vorticity maps of a General Electric J85 undergoing stall 
due to pressure inlet distortions. It was shown in Section 
IV that single gas, incompressible flows with temperature 
distortions can be treated as pressure distortion problems. 


thus there is little loss of generality in using Farmer's 
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a. work. Figure 13 shows a circumferential vorticity contour 
map computed by Farmer indicating the radial of interest. 
Since the J85 engine has annular cascades and rotors 

with varying stagger angle and tip clearances, three addi- 

tional and very coarse approximations were made. First, 

the annular cascade was approximated by a rectangular cascade. 

Second, the varying stagger angle was replaced by a mean 

stagger angle. Third, blade tip losses were excluded. 

These assumptions were necessary to remain within the con- 

fines of the Squire and Winter expression. 

. The solution of Poisson's equation was obtained by using 
numerical methods on a digital computer. Appendix G presents 
_the problem formulation and method of solution. There are 
presently available several computer routines which can be 

: employed in more general flow situations. Hawthorne and 

Novak [Ref. 42] present a brief summary of the routines 


available at the writing of their article. 


D. DISCUSSION OF RESULTS 
The stream function (y) distribution is listed in Table 
Ii. Examination of the change between the first and second 


column is equivalent to considering the spanwise induced 


velocity w= - a - Similarly, the change between rows 
; is equivalent to considering the induced downwash v = , 


Qualitative analysis of these data results in a streamwise 
vorticity distribution as shown in Fig. 14. 
These vorticies have two effects. First, the Biot- 


Savart law requires that their influence be felt by the 
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Figure 13.° Circumferential Vorticity Contours, J85 


(adapted from Ref. 13) 
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THE PSI DISTRIBUTION 


(READ ROW-WISE. FIRST VALUE IS FOR Z= 0.111 AND Y= 0.111) 
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Figure 14. Cascade Secondary Vortices 


Due to W4 Inlet Distortion 
onl 





circulation round the blades of the inducing passage. This 
is shown in Fig. 15 where their induced downwash is plotted. 
Note the decreased downwash over the outer portion of the 
blade. This characteristic of the flow decreases that sec- 
tion's angle of attack. The second effect occurs as the shed 
vorticies are swept downstream where they will influence 


the succeeding stages. 





Figure 15. Induced Downwash 
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VII. SUMMARY 


The phenomenon of TFG distortion has been discussed from 
the vorticity viewpoint. This approach provides insight 
into the aerodynamic influence associated with TFG distor- 
tion. A pure temperature distortion problem was analysed 
using a temperature map recorded just prior to engine stall. 
An order of magnitude evaluation of the vorticity equations 
reduced them to their essential terms. feaitueeien of these 
terms provided data from which vorticity maps were construc- 
ted. It was found that by themselves the maps provided some 
pertinent results which enabled limited analysis to be made. 
However, to derive full benefit from vorticity analysis, 
time history and engine performance data were found necessary. 


The technique used to develop the vorticity .maps was 


strictly manual in nature. There seems to be great promise 


if this process could be generalized into a computer program. 
The main benefit would be labor savings which in turn would 
allow time history temperature or pressure maps to be con- 
verted into time history vorticity maps. 

The discussion surfaced various voids of knowledge. 
For example the effects of unsteady flow and cascade second- 
ary flow in a TFG environment are seemingly unknown. The 


equations of vorticity require flow density or concentration 


a3 








maps which are generally not being considered in current 
TFG studies. 

The idea of the aerodynamic transfer function summarizes 
one suggested use of the vorticity approach. This is be- 
cause if one combines the compressor stage's design charac- 
teristics with its aerodynamic response characteristics, 
the stage transfer function is defined. Vorticity is the 
initial input. Figure 16 illustrates this point. Do 
is the vorticity impressed on the engine inlet. The stage 
transfer function is G. As a result Do_y is transformed 


into D ~ the distortion from the first stage being 


1-2 
impressed on the second stage. Note the feedback loop which 





Figure 16. Block Diagram for 


Aerodynamic Transfer Function 


tit is recognized that collection of flow density data 
presents a problem. In most experiments this problem has 
been glossed over by considering total compressor or engine 
response. The fact is this does not attack the problem of 
determining what is aerodynamically happening to the blades 
or cascades. 
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. would likely be present in subsonic flow. Granted, this is 
conceptually easy to visualize and extremely difficult to 
implement. But with vorticity one has an aerodynamically 
compatible input with which to visualize and analyze inlet 


TFG distortions. 


De 





VIII. CONCLUSIONS 


As a result of this study, the heretofore abstract inlet 
temperature map is now a versatile tool in the analysis of 
temperature-induced instabilities in an axial flow compres- 
sor. The vorticity concept provides answers to questions 
like "How is the blade circulation being altered?" and 
"What sort of secondary flow effects can be expected?" 
Answers to both questions can be found through respectively 
radial and circumferential vorticity. Mathematically a 
third component,.namely axial vorticity, ean be obtained 
from a two-dimensional map. 

To optimize this technique, additional efforts are re- 
quired. The vorticity equations should be generalized to 
cope with flow situations involving chemical reactions or 
phase changes. The effects of TFG on current unsteady 
airfoil theory and cascade secondary flow eneeiey require 
investigation. Flow density mapping procedures need to be 
developed. The latter two areas would not only benefit 
the vorticity concept but would enhance the general under- 
standing of the TFG problem. 

The vorticity concept leads to the idea of an aerodynam- 
ic transfer function. In turn, this transfer function is 


@ move in the direction of the ultimate goal ~ prediction of 


flow of an arbitrary compressor with arbitrary inlet flow. 





APPENDIX A 


Derivation of TFG Vorticity Equation 





Figure A-l. Natural Coordinate System 


(adapted from Ref. 45) 


Figure A-l is an example of the natural coordinate system, 
This system will be used to develop a general vorticity 
relationship involving pressure, temperature, and foreign 
gas. The development will follow that of Liepmann and Roshko 
[Ref. 45] with the modification of presence of foreign gas. 

When a steady, two-dimensional, inviscid,adiabatic flow 


of a non-homogeneous mixture of perfect gases is considered, 


the equations of motion are 

: eUAn = constant, (A-1) 
Pe. gp oll e 
35 eu ss es 


Of 








and 
oP _ US . _ 
oT PR - 
By definition 


h, = h + 


2 
‘ U (A-4) 


nie 


Equation A-l represents continuity for a unit width 
streamtube. Equations A-2 and A-3 are respectively s and n 
He mew Putt balance. Equation A-4 is the Nenergy" equation 
where by definition stagnation enthalpy is constant along 
any streamline. 

th 


af ms represents the mass per unit volume of the i 


perfect gas at any given point,then let 


n 
M= } mM,» (total mass/unit volume) 
1 
x, = m,/M ; (mass fraction) 
g meee 
i 9s? 
and 
f = 
i on 


The first two equations are definitions. The second two 
equations are arrived at by assuming gradients in the n and 
s direction of the mass fraction of each gas. 

For mixtures of perfect gases, the second law of thermo- 


dynamics holds that 


7 

S$ = qT as 

5 =) om, (ey Inge Ry the.) (A-5) 
1 i r r at 





where: 8 = entropy of mass M 


s = specific entropy 


Subscript r = 


reference condition 


Superscript bar used to delineate between entropy (8s) 


and the natural coordinate (s) 


In terms of specific entropy, equation A-5 can be rewritten 


a6 
— n — 
5 = ) X85 
1 
Thus 
95 n oS. _, Sm, 
on } Xs On Si; on 


By oT oP, 
Y [xy |-g> ga 7 PL On (afd 
1 
(A-6) 


In most practical cases involving distortion the number 


of perfect gases can be considered to be two - air and the 


foreign gas. 


- relationships apply 


X, +X, = 1 
P, +P, =P 


V,=V,=V 


Considering two species of gas, the following 


It is emphasized that chemical reactions or phase changes 


are not included in this formulation. 


Equation A-6 can be expanded and rearranged so that 


De 








_ 1 
= <= = ob — 
3a f(s, s,) + a (x a, Co. 
x R_ oP x R 9p 
11 1 2 2 2 
as: on P on 
By oéfinition 
x,c XC = C of the mixture and 
Pp, Po 
x,1R, Halts 1 
Po . ev 
which leads to 
os _ ~ = oT 1 aP 
acral TR, Cs, iS) Cy I~ Pp on 


Similarly, it can be shown Eq. 


as 


ae = 1 
oe = + eee acc) ae 
7 Tg (8, s,) C 5 OS 


Ss 


For a chemically inert, perfect gas mixture where h 


U 


t 


then 
o af = ah 
p on n 
From Eq. (A-4) 
Os ee. 
on dn 
Thus 
: ar dh 
HD om 7 dn 


By introducing Eqs 





(A-7) becomes 


p os 


aC 
Gf tae 
on 
au 
on 
oC 
pee ae 
on on 


(A-~3) and (A-10) into Eq. 


dh, 3C 
ME ae ae = U 


(A~8) 


(A-10) 


In natural coordinates the flow vorticity is defined as 


ye BW _Y 
~ on R 
leaving 
oS Ses 8 ass 
T an = Tf, (s, = s,) + an = 7 an - Uw : (A-11) 


The same procedure applied to Eq. (A-9) results in 


a5 aC 
7 ae Tg (s - s.) - 7 5 


Note that dh, /9s = 0. Assuming no diffusion, condensation, 


or chemical reactions 


and 


Is _ 
te ee ee ; (A123) 


Equation (A-12), which is the result of steady, inviscid, 
adiabatic flow, indicates that entropy along a streamline is 
a constant. 

By subtracting Eq. (A-8) from Eq. (A-11) 


dh aC 


_t mm _P_ or _ 1 oP 
dn . an Cy on se p on 
This can be reduced to 
dh 
oe ge ee Bh 
Uy aa + > On 5 (A=12) 


which is the desired result. 
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APPENDIX B 


Vorticity in Incompressible Flow 


Consider the system shown in Fig. B+1l. Assuming the flow 
is inviscid, adiabatic and incompressible, Eq. 5 can be 


applied. The next step is to write the corresponding 


t t 
Pr © 
—y, ap) 


Figure B-1. Flow Model 





scalar equations. In this case the cylindrical form was 
chosen. A suitable scheme should be employed to.nondimen- 


sionalize the terms. The result might look like 


oP 2 dU 
oi ov. Uo sg 
Upe, 7 Un p or 2 sr Ot : (B-1) 
oP 2 oU 
7 i ee BE BIS es 8 - 
U4, — Uw, = 6 36 r36.0CU Csi‘ (B-2) 
and 
oP 3oU 
o Oe . Ue elie Zz 
UW, Vets ~ Fae 3 a2) ‘ (B-3) 


where U is the non-dimensional velocity, w is the nondimen- 


sional vorticity and the subscripts indicate the component 
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direction. From the initial assumptions 


9P | 
—H=9 anf “88 2'p 
a2 OZ 


One would have a total pressure map (or series of maps) 
representing the inlet conditions at some event of interest 
i.e. stall. From this map a maximum and reference pressure 


(P,) could be selected to define an order of magnitude e 


It becomes a matter of using the available data and/or 
analyzing the physical flow situation to arrive at orders 

of magnitude for the components of the equations. For 
example,if there is small duct swirl (w) and small acceler- 


ation in swirl then one might assign the following estimates: 


aU 
lL. w.=0(e); 2. Uj=O(e); and 3. spe =0(e). ° From physical 


considerations the radial component of the flow should 


.. generally be neglhigible with small acceleration. Thus 





U,,=0(e?) and aU,,/at is O(69. The axial velocity U,)can be 
assigned O0(e€°). The gradients in radial and circumferential 
pressure can be determined from the pressure map. This 
procedure would be carried through until all terms had estimates 
as to their relative order of magnitude. By neglecting terms 
and products of terms of order € or less, the equations might 


look like 











UM, = or 060 or 26 : 
and 
ou, 
gee = 


Conceptually, these equations now express vorticity in terms 
of the flow variables. A more detailed procedure is carried 


through in Ref. 13. 
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APPENDIX C 


Order of Magnitude and Sign Verification of Eq. (6) 


Consider the model shown in Fig. C-1 (a). 





(a) (b) 


Figure C-1. Flow Model 


Assume: (a) Steady, one dimensional, insentropic flow from 


1+ 2. 
= 5 ° ot © = 
(b) T. 520° R T. 620° R Vs Lad 
1 1 
YA = P. = 1 atm Pa = P, v= dee 
1 1 2 2 
= ft - lbdf ; 
R= ; Ea ina ae = 
a 2343 lbm - °R “a oe? 
e Seoee. Fee 
SS ? lbm - °R 


(ec) Linear variations in temperature and mass frac- 


tion as depicted in Fig. C-1 (b). 








From isentropic tables for a 1.4 and M, = 0.5 
; 2 
P T 
oe 84 a5 
= 0,943 = 0.952 
P. qT, 
a a 
2 2 
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For isentropic flow, Eta, = Pta, = 1 atm. and "es, = Tea, 


0.843 atm. and Ta, = 495° R, 


=T, = 520° R. Therefore, Pa, 
1 


From isentropic tables for y, = 1.3 and Ps /?ts = 0.843 atm. 
2 


M = 0.52 


i 





0.961 . 


Therefore, Ty, (.961)(620) = 595° R. For speed of sound, 


a; 
a (ye, RT)*. Therefore 


a, = Wt. A200) (seo) 49D 1090 ft/s 


a 


1460 ft/s 


tl 
oN 
bw 
Ww 
ty 
ine) 
NM 
Oa 
On! 
XO 
WI 

lt 


a 
iS) 


The velocities are 


Cc 
1" 
S 
© 
nt 


(0.5)(1090) = 545 1t/s 


and 


(om 
tt 
=S 
a) 
tl 


(0.52)(1460) = 760 ft/s 


Since the model is for steady planar flow with constant 


static pressure, Eq. (6) reduces to 


dh, ah 


~wU = dy dy 


Recall that Us = X Cp tk cp, and f, = dx,/dy. Rearranging 


the above in terms of specific heat and temperature leads to 


aT t 
JV = f (T - T Cc - Cc + See See 
o a 4) 6 p ee cal 


in linear and finite difference form 
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7... #2, 
UU} {U4 x -0)| 7. * t, ts 
Ay | AY P, Ps 
Ty =T 
T —T t oe ec +e 
+ a, Ss, = a, o5 Pa Pa 
Ay AY = 
which reduces to 
(U2 = U2) = 2c. (Tt. = 7 ) + 2c (T a. 4) 
2 - Pa ae “a, Ps oS S) 


Evaluation with values previously calculated 
2,8 « 10° ft*/sec* = 3.0 * 10° ft*/sec* 


The sign is consistent, and the error in this crude approxi- 


mation iS 6.5%. 
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APPENDIX D 


Analysis of Products of Combustion 


On a mass basis the chemical equation for the combustion 


of hydrogen in theoretical air is 
2.02 lbm Hy + 16.0 1lbm 05 + 106 lbm No 


# 18.02 ibm H,0 + 106 lbm Ny eo (D-1) 


The N, term accounts for nitrogen as well as the other inert 
gases in air. From the data in Ref. 39 the mass flow rate 


of H, for a 150° R change across the heater was computed to 


2 
be 0.27 lbm/s. The rated air flow of the engine was 186 


lbm/s [Ref. 39]. If air into the heater is assumed to be 
dry, then, based on a one second interval, Eq. (D-1) modified 


for excess air becomes 
0.21 lbm H, + 39 lbm 0, + 147 lbm No 


= 243 lbm HO + 36.80 1lbm Dy, + 147 Ibm N,. (D-2) 


From the results of Eq. (D-2) the specific humidity is found . 
to Be 1.3%. Based on that small value, two conclusions were 
reached: 1. The effects of this humidity on engine perfor- 
mance are negligible [Ref. 46]; and 2. The inlet flow 


can be considered to be free of foreign gas (water vapor) 


gradients. 





APPENDIX E 


i Heat Addition in Subsonic Flow 


The following equations were developed by Tsien and 


Beilock [Ref. 47]. For a point heat source in subsonic flow 


2 See eee aa 
U omy PB x2+p2 y? ‘ e ) 
é (¥-1)Q y 
vi = es Ue cere ssearnic Damen a 
ony PB xiapiy? (E-2) 
mn and 
(y-1)MQ x 
% = SRE aE e El - 
P enap x2+p2y? : (E-3) 


In this notation x and y are analogouS to z and r, respec- 
tively. Here * indicates a change to the flow variable 
between the heat source (point 0) and some point (point 1) 
down (or up) stream. For example, U, 0) Um 5 a +P*® etc. 
Rise g=1-M?, a is the speed of sound and Q is the line 


source heat addition (BTU/ft-s). If one imposes the require- 





ment that P*=0, then for the non-trivial case x must be 
large. So, when P*=0, V*=0 and U*=0. Physical considera- 
tions lead to the conclusion that Ug and Uy, are of the same 
order as V*. Thus they will necessarily be small. 

For the aU L/at! and 9U'/gt' analysis, Eq. (E-2) can be 


differentiated with respect to time 


aig | Fat #8] 


a5 es 
at 2m 2yg2y2 L T Bt t 
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Again, for large x, it can be shown that av¥/at=0. Using 
the same rational as applied to U4 and UL aU Tot! and 


aU, / ot are small. 
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* APPENDIX F 


Order of Magnitude Analysis for a Heat Ramp 


The system shown in Fig. F-1 is a constant area duct with 
heat addition at station 1. Assume the heat source is capa- 


ble of producing a temperature ramp of 7750° R/s. Suppose 


t 
7 


® 





Figure F-1. Heat Ramp Flow Model 


one had a case of Rayleign flow with an impulse change in T, 
: 1) 
and a corresponding impulse change in P. as shown in Fig. 


F-2(a) and (5). 


a | oe een 


time time 
(a) (b) 


Figure F-2. Impulse Changes in the Flow Variables 


ti 


The following assumptions are made for the initial con- 
ditions: 1. P=1 atm; 2. M=0.2; and 3. U = 200 ft/s. 
If the flow particles are convected at a velocity U, 
then changes in pressure propagate at Uta. Changes in tem- 

perature due to particle convection propagate at U. 
From the assumed initial conditions the following addi- 


tional initial values can be computed 


a= U/M = 1000 ft/s , 
= 27 = Hig 
T = Yer 416° R . 
and 
= Yo aes 2 ° 
T, 1, (i + 5 M*) H15° F. 


At the assumed Mach number the change in Tt. can be approxi- 
mated by the change in T. Thus, between t = 0 and t = 0.01 
seconds, AT, = 77.5° R, etc. Computation of changes in the 
remaining flow variables during the first 0.01 seconds can 
be accomplished through the use of thermodynamic influence 


coefficients 


It 
Wt 
oO 
ine) 
o 
un 
SG 
" 


= 240 ft/s 


dP 2 dT 
ape = - UES —+ = - 0.005 or P, = 0.99 atm . 
t lt S=M’ T 


In a similar fashion, progressive changes can be computed. 


Results for the first 0.02 seconds are listed in Table F-I. 


Te 


‘This is an approximate analysis; the correct method would 


involve use of wave techniques outlined by Rudinger [Ref. 48]. 


Table F-I. 


Time Changes in Flow Variables 





£2 1000 416 19 1.0 
,22 1090 »=—- 4g 4g7 0.99 
.24 1172 572 5T5 0.98 





It is observed that with respect to tame the ehanges in 
velocity and temperature are in order of magnitude ten times 
that of pressure. 

From the values in Table F-I, Fig. F-3 was constructed. 
The slope of the lines are the reciprocal of the transport 
velocity. In this case the changes are with respect to dis- 
tance. Again it is obvious that the changes in pressure are 
negligible. Hence, it can be concluded that for a tempera- 
ture ramp process the resulting gradient in pressure is 


negligible. 


13 


lines of constant temperature 


lines of constant pressure 


Time 


(s) 





1 2 3 4 


Distance - 
Care 


Figure F-3. Transport Characteristics 
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APPENDIX G 


Numerical Solution of Poisson's Equation 


Poisson's equation is an equilibrium problem involving 
an elliptic operator. The boundary value problem is formu- 


lated by a governing equation 


where V* is the Laplacian operator, y is the stream function, 
e is the turning angle of the flow and oa is the inlet veloc- 
ity distortion. The boundary conditions are p = 0 at y = 0, 


z= 0 and2Zz=2z . In using inlet circumferen~ 


Y= Vmax ? max 


tial vorticity the governing equation is modified to 


vey =~ 260, 


since by definition We = + a 5 

The solution technique involves "discretizing" the domain 
and writing finite difference approximations for the con- 
tinuous derivatives. These are the standard techniques de- 
scribed in most numerical methods texts. 


In matrix notation the general result of the finite dif- 


ference operations is 
[D]{v} = - 2e{ug} . 


where [D] is a matrix of finite difference coefficients. The 


solution is obtained by the linear algebra notation 


{y} = - 2e[DI""{ug} 


He 





A unique feature of [D] is that it is a banded matrix. 
This feature can be exploited thru the use of banded ma-— 
trix solving subroutines. One such subroutine used in this 
problem is GELB which is found in the IBM Scientific Sub- 
routine Package [Ref. 49]. 

To test this program, Squire and Winter's [Ref. 37] ex- 
perimental data was used. The results are shown in Fig. 


G-l. 





w = velocity in z direction. 


U. = mean velocity of flow approaching bend. 


abscissa = inches from inner vane. 





Squire and Winter's measurements. 


--O-- Numerical solution. 


Figure G-l. Comparison of Squire and Winter's Experiment 


with Numerical Solution 


* The domain of Farmer's [Ref. 13] J85 problem was modified 
from Fig. G~2(a) to Fig. G-2(b). The stagger angle which 


varied from 22° 25' at the hub to 52° 59' at the blade tip 
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| B50" 


Babe 


Lek? -| 
(a) (b) 
Figure G-2. Actual and Modified Domain 


for J85 Solution 


was replaced by a mean stagger angle of 37° 45°. The turning 
angle e« was computed by Farmer to be 19°. Values for 


were computed using the equation obtained from Farmer's 


paper 
-2854., f 
ee 5 
P'/9.i7P' — 8.16 or! 


Values for 3P'/3r! and P' were obtained from Fig. 13 along 
radial I. To obtain the component of Wg parallel to the 
inlet plane (Fig. G-3), the values from Eq. (G-1) were multi- 
plied by cos(37° 45'). The incremental step size was chosen 
to be 0.111 inches which resulted in a matrix [D] of dimen- 
sion 49 x 9. 

The program appears after the appendices. Since GELB 


is a standard IBM subroutine it was not included. 


(7 








° ; Figure G-3. Components of w 


8 
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ANAMHAAANRDAAAAAAANANAAANAAON 


aAaAM AAD AAO AON 


AAD AHO 


QAM 


os 


LB) »,DISP=SHR 


SCI a eRe aK a akc Hi OK AC aK OE OK IK HEE Se Se oe ie feo ae a OK aK ak a a ae ae ee 2 kc a Be a ae ea ae kee 9K 


NUMERICAL SOLUTION OF POISSON'S EQUATION. 


DESCRIPTION OF PARAMETERS 
BD —- VECTOR OF DIMENSION MA CONTAINING COEFFICIENTS 
FROM FINITE DIFFERENCE SCHEME. 
FUNC - VECTOR OF DIMENSION M CONTAINING ~2*EPS*WTH. 
WORK ~ VECTOR OF DIMENSION MUD USED AS A WORKING 
WIH - VECTOR OF DIMENSION M/MUD. 
M = L/STEP 
STEP - INCREMENT SIZE. 
EPS - FLOW TURNING ANGLE. 
EPSL — A TOLERENCE FOR GELB. 
REMARKS . 
1. FOR DESCRIPTIONS OF THE PARAMETERS MUD, MC MAy 
ME AND N SEE GELB. 
2. ROW IS ASSOCIATED WITH HEIGHT Z. 
COLUMN 1S ASSOCIATED WITH WIDTH Y. 
ake ae ai ake ois ake ok ei aK ie a aK he se ok sake ak te ah ak i oe ai ok Re OK AE Ne ke fe ofc of ee RC akc OK ae aC aK OK ac aK a aK 2K as 2K a Co 
DIMENSION 0(8334) ,FUNC (441) pWORK(19) »WTH(49) 
DATA MzMUDsMCy MAy MEyN/ 4419 99197 8334, 82899 1/ 
DATA STEP EPS EPSL/LeL1OE—Ol 91 .0E-06 4363160E-01/ 


KZ=M/MUD 

READ VALUES FOR OMEGA-THETA. 
READ (5,100) (WTH(J)sJ=19KZ) 
ECHO CHECK. 


WRITE (6,101) 
WRITE (6, 107} (VITHC IS) J=12KZ) 


FORM THE FUNCTIONAL VECTOR ‘FUNC’. 


Oc I 
li 


TMNQCReOOe 


—2 .O*EPSL¥WTH( JJ) *COS( -6555) *STEP* STEP 
' OF EACH ROW. 


EP : 
UNC{J) yJ=1 M9) 


el 
DPRAD A CHeCHTN 
eit eet tet oY) | Coa 


SEs 


FORM THE COEFFICIENT VECTOR ‘D°. 
CALL LOAD (WORKsDsMsMC »MUD yME_MA) 
SOLVE FOR THE STREAM FUNCTION. 


CALL GELB (FUNC,DsMyN,MUD,»MUD,EPS,TER) 
DF A ITERLEQ<-1) CO 1G 39 


PRINT STREAM FUNCTIGN. 
WRITE (6,105) STEP,STEP 


1 
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LOAD LOWER VEC WITH UPPER MOLV- 


ZERO REMAINING MA — ME SPACES OF VEC. 


INITIALIZE SIDE BOUNDRY CONTITIONS TO ZERO. 


~éMO>Y~>oxwW 


ina) 
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